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A general expression is derived for the dipolar NMR spin-lattice relaxation rate 1/Tj of a system exhibiting Brownian 
dynamics in a discrete and finite configuration space. It is shown that this approach can be particularly useful to model 
the proton relaxation rate in molecular rotors. 


I. INTRODUCTION 

Artificial molecular rotors are attracting much interest 
as a possible route to achieve a controlled motion at the 
nanoscale 1 S A lot of progress has recently been made in 
solid-phase based systems as molecular engineering permits 
to conceive crystalline environments where such rotors can 
move more-less freeljEEl For molecular-sized rotors, even 
though a preferred direction of motion can be induced by ap¬ 
plication of an external driving field, the main component of 
motion generally consists of random reorientation due to ther¬ 
mal vibrations of the crystal lattice. 

The experimental studies of such systems need tools sensi¬ 
tive to the dynamics of motion as complementary to average 
structure determination by X-ray diffraction. The latter gives 
mean positions of the atoms and the amplitudes of their av¬ 
erage displacement from the mean position, however it does 
not give information about the timescale of the displacement. 
Though, from the experimental point of view, it is the study of 
the dynamics of the thermal motion which provides the main 
source of information about the interaction between a rotor 
and its environment. Among various experimental techniques 
that have been used for this purpose 1 ’“NMR is probably the 
most powerful one as it can probe the dynamics with atomic 
resolution. In crystalline-environment solid-state NMR tech¬ 
niques can be applied similar to standard liquid-state NMR 
methods that are commonly used to study dynamical pro¬ 
cesses in solution. 

For crystalline systems most often basic NMR techniques 
are used, studying either the spectral shape or Tj relaxation. 
These two approaches differ by the range of frequencies of 
the motion that are probed. The first relies on line-narrowing 
effects appearing when the frequency of the motion is com¬ 
parable to the inhomogeneous NMR linewidth resulting from 
the local field anisotropy, whereas the spin-lattice relaxation is 
sensitive to motion at the timescale of the Larmor frequency. 
For example, deuterium NMR provides a very good probe as 
the lineshape is affected by strongly anisotropic quadrupolar 
interaction and simulation of 2 H NMR spectra has been of¬ 
ten used in deuterated rotorJSB 15 . However this technique 
needs deuterating the sample (proton linewidth in solids is 
dominated by the dipolar proton-proton interaction which is 
less affected by a constrained anisotropic motion and the line- 
shape is not so easy to simulate). For protonated samples 
: H Ti relaxation provides instead a simple tool to study the 
dynamic ^ 1 ° I 13 I 16 1 17 [ 

Proton spin-lattice relaxation by thermally activated mo¬ 


tion, where the relaxation process is induced by a random 
modulation of dipolar interaction, is usually modelled using 
the well known Kubo-Tomita (KT) theory-®. However the 
standard KT formula was derived for the simplest case of 
isotropic reorientation with dynamics described by a single 
correlation time, and for more complex systems the problem 
is to know the distribution of the correlation times and their 
relative weights in the 7j relaxation rate. Such models were 
developed for polymers and biological materials 19 20 . This pa¬ 
per deals with relatively small systems that can be character¬ 
ized by a finite number of correlation times, such as molecular 
rotors. 

In the absence of an external driving field the thermal mo¬ 
tion is determined by the torsional potential. The equilibrium 
positions correspond to the potential wells and the nature of 
the motion depends on the height of the barriers A between 
neighboring wells with respect to the thermal energy /cl®. If 
the latter dominates the kinetics of the motion is a Brown¬ 
ian random walk that can be modeled as a diffusion process 
(i.e. the probability of finding the rotor at a given position in 
function of time is governed by diffusion equation) where the 
diffusion constant depends on the effective friction and the in¬ 
ertia of the rotor. The model presented in this work present 
deals with the opposite case A >> kT, probably more of¬ 
ten encountered in crystalline environments. Here the motion 
is hindered by the barrier, in this condition the motion mainly 
consists of thermal librations within a well of the torsional po¬ 
tential but the rotor also occasionally undergoes thermally ac¬ 
tivated hops between adjacent wells. The thermally activated 
hopping rate T is usually modeled by the Arhenius law 

r = woexp— A/kT (1) 

where the attempt frequency ojq is close to the librational 
frequency 121 . The condition A >> kT implies T << wo 
therefore, for small rotors in crystalline environment the spec¬ 
tral distribution of the characteristic frequencies of the mo¬ 
tion is concentrated in two distinct regions. This is in con¬ 
trast to disordered polymeric materials or biological macro¬ 
molecules where the motion exhibits a rather continuous fre¬ 
quency spectrum 19 20 . Here the libration frequencies are typi¬ 
cally of order of 10 13 Hz which is several orders of magnitude 
higher than the typical NMR frequencies therefore NMR is 
rather insensitive to this spectral component. On the other 
hand, the hopping rate V can match the nuclear frequency at 
sufficiently low temperature leading to an efficient Xj relax¬ 
ation process. 

The fact that NMR is only sensitive to the slow spectral 
components of the motion is an interesting feature since it 
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means that, as far as NMR is concerned, only the positions of 
the potential minima and the height of the barriers matter but 
the exact form of the potential is not relevant. This permits 
to make very simple and robust models for NMR T\ relax¬ 
ation. The first such model has been proposed by Bloember- 
gen et al. 22 assuming that the random fluctuations of the local 
field due to hopping is described by a simple correlation func¬ 
tion of the form exp —Ft. The Kubo-Tomita formula-^results 
from the application of this model to the case of relaxation by 
random modulation of the homonuclear dipolar coupling (see 
next section). 


The simple exponential form of the correlation function 
is justified for hopping through a single potential barrier be¬ 
tween two equilibrium states. This theory applies well for e.g. 
hopping between two conformational states of a molecule. 
The purpose of this paper is to complete this theory with corre¬ 
lation functions calculated for an arbitrary number of degrees 
of freedom, wells and barriers. Examples that will be studied 
are rotors having different non-equivalent positions and a pair 
of coupled rotors. Inducing and detecting correlated motion 
in systems of coupled rotors is also one of the current chal¬ 
lenges in the field, the studied example may provide the back¬ 
ground for understanding the results of a recent work where it 
was found that the 7) relaxation in a pair of coupled rotors is 
dominated by two different activated processes 17 23 . 

In the literature devoted to molecular rotors, the term 
’’Brownian motion” is often reserved to the case of purely 
thermal and unhindered motion where the thermal energy kT 
is dominant, that is the case A << kT. But of course even 
in the opposite case considered here the hopping through the 
potential barriers is also a Brownian process. In this regime, 
once we forget the librational degrees of freedom almost invis¬ 
ible to NMR, we are left with a system contained in a discrete 
configuration space (that is where space coordinates only take 
discrete values) defined by the set of wells, exhibiting a Brow¬ 
nian dynamics. 

The theory of the Brownian diffusion for rotating objects 
and the calculation of the correlation functions involved in 
NMR relaxation can be found in a number of works. For ex¬ 
ample, the textbook case of fully isotropic rotation (diffusion 
on a sphere) is treated in A brag air. 23 , and a more general case 
of anisotropic rotation has been considered in 25 29 . The calcu¬ 
lation presented here follows the same philosophy except that 
it uses somewhat different mathematical tools adapted to the 
finite discrete space and will therefore be detailed for com¬ 
pleteness. Nevertheless it should be stressed that the theory 
presented here is in some sense more general because it is not 
merely restricted to a simple geometry of motion such as rota¬ 
tion. The discrete configuration space has no dimensionality 
in the usual sense (its topology is entirely defined by the set 
of hopping rates between different states), therefore the same 
formalism can be applied for any system that has a finite num¬ 
ber of states, such as a pair of coupled rotors. 

The paper is organized as follows. In the next section we 
derive an expression relying the correlation functions to the 
hopping probabilities and apply it to calculate the dipolar re¬ 
laxation rate 1/Tj. In section III we present two situations 


where this approach permits to go beyond the Kubo-Tomita 


analysis and calculate the variation of the relaxation rate 1/Tj 
from microscopic parameters. 


II. THEORY 

A. Basic formalism for motion-induced T1 

In this section we recall the basic formulas for the relax¬ 
ation rate of a pair of nuclei induced by random modulation 
of their mutual dipolar interaction. After A brag am 2 -, the gen¬ 
eral form of interaction hamiltonian can be written as 

'Mint (<) = £ F^\t)A^ q) ( 2 ) 

Q 

where F^(t) = F^*{t) are some random functions of 
time related to the positions of the nuclei and A <q> are spin 
variables. The index q is usually related to the change in the 
z component of the total spin for a given relaxation process. 
Note that other mechanisms of relaxation induced by random 
motion, such as modulation of the quadrupolar interaction or 
the shielding tensor anisotropy, can be cast into this general 
forrrPS 27 29 . Here we focus on the relaxation by dipoledipole 
interactions between identical spins, which is often the domi¬ 
nant source of relaxation in molecular rotor systems carrying 
protons (for the other cases we need to replace the Eqs.Q 
below by the appropriate expressions). 

For this interaction, considering two identical spins 1\ and 
I 2 we have 


= -\ 1 2 h{h z h± + h±h z ) 

A (±2) = — | 7 2 W 1± /2± 

F (±1) = isinflcosfle^ 

F (±2) = A sin 2 fle* 2 ^ (3) 

Here the vector r = (r, 9 , (j>) defines the relative positions of 
the spins and the time dependence of F <JI) (t) comes from the 
random modulation of r(f). 

Then the relaxation rate 1/Tj of the pair of spins is related 
to the spectral density of fluctuations of F ,q> ( t ) at the Larmor 
frequency uin- Let us define the correlation function 

$,(t) = (f’ (9) (r)F (9) *(t-fr)) T (4) 

then the spectral density is 



and the relaxation rate of the spin pair is given by 

1/Ti = (7 (Ji(wjv) + J 2 ( 2 u N )) (6) 

with C = | 7 2 ?i/(/ + 1). For an isotropic random ro¬ 
tation of the vector r characterized by a single correlation 
time t the spectral densities ar^ 2 - 824 Ji(cc) oc 1+ F . T 2 and 
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FIG. 1. The discrete model is defined by a set of potential wells E l 
and barriers Ay between them. 


h{u) oc 1+ % T $ and the Eq. § gives the well known Kubo- 
Tomita formula. 

Note that the description of the system based on the hamil- 
tonian Q works under the assumption that the motion can 
be treated classically and that the orbital and spin variables 
are separated. This does not hold in the case where the 
system may undergoe quantum tunneling between configura¬ 
tions having different spin states. A well known example is 
the quantum tunneling observed in CH 3 groups, this process 
yields, at low temperatures, an additional relaxation channel 
that is not captured within the hamiltonian (| 2 ]PSI. 


B. System dynamics: diffusion in a discrete space 

At present we consider a system defined by a conforma¬ 
tional potential exhibiting some number of wells of energy 
Ei separated by barriers Aas sketched on FigJT] The hor¬ 
izontal axis in FigJT] represents schematically the degrees of 
freedom that may have more than one dimension (e.g. two 
angles for a pair of rotors, as will be considered in the next 
section). As stated in the introduction we consider the case 
A ij » kT, that is when NMR mainly senses the slow spec¬ 
tral component of the motion due to thermal hopping between 
the wells. Under such conditions the dynamics relevant for 
the Ti relaxation can be modelled by diffusion in a discrete 
space defined by the set of states {Ei}. Let rii be the popula¬ 
tions of these states, it will be convenient to use the convention 
rii = 1 i.e. population is defined as the occupation proba¬ 
bility. We denote as the populations at thermal equilibrium 
(then rii ~ exp {—Ei/kT)) and as T'y the probability of hop¬ 
ping i —> j per unit time (i.e. the average hopping rate), for 
a thermal process, according to 0 this probability is propor¬ 
tional to exp(—A ij/kT) where the barrier Ay is defined in 
FigJT] Since Ay — A,.,; = Ej — Ei we therefore also have 
T i:i = Tji exp {-{Ej - Ei)/kT ). 

The time evolution of the system from a non-equilibrium 


state (relaxation) is governed by the master equation for pop¬ 
ulations: 


drii 

dt 


'y [ (j) Tijn^t) 


(7) 


This equation can be rewritten as 


dni v—■v . 

-to =-^2 A a n i( t ) 

3 

with 

4 = S~ T Ji 

jE t r ifc for i = j 


Then, in matrix notation with 11 = A 

master equation is 


( 8 ) 


(9) 

{Ay} the 


dn 

dt 


= — An 


and has a general solution 

n(f) = e" A 'n(0) (10) 


Therefore, the eigenvectors of A define the ’’diffusion modes” 
having different relaxation times related to the corresponding 
eigenvalues. 

As an example, consider a rotor having N equivalent po¬ 
sitions (i.e. all states having the same energy and separated 
by the same barrier) and a uniform probability T of hopping 
between adjacent states (i —i + 1 and i —> i — 1). Then the 
master equation reads 

= T {n i+1 + rii-i - 2rii) (11) 

dt 

Note that the right hand side of this equation can be written as 
T ((n .; + 1 — rii) — ( rii — rii- 1 )) showing that it is obviously a 
discrete equivalent of the one-dimensional diffusion equation 

dn{x,t) d 2 n{x,t) 

dt dx 2 

For example, for N = 3 (equilibrium positions at 0, 120 and 
240 degree) the matrix A is 



The corresponding eigenvectors are: n\ + 77.2 + ^3 with the 
eigenvalue 0 (the zero eigenvalue is always present because 
the total population number is conserved: n, =const), and 

n-i +ri 2 — 2ri3 plus permutations with the eigenvalue 3T. Thus 
1 /3T is the characteristic relaxation time of this system. 

For a general N it easy to show that the equation (111 
transforms to diagonal form with eigenvectors of the form 
rii ~ cos ki yielding the eigenvalues 2T(1 — cos k) where 
the allowed values of the wave vector k are determined by pe¬ 
riodic boundary conditions cos (ki) = cos(fc(z + N)), e.g. for 
N = 3 these values are k = 0, 27r/3 and 47r/3 corresponding 
to the eigenvalues 0,3T,3T. 
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Calculation of the correlation functions. Let F be any 
function of state of the system, for a discrete configuration 
space F is defined by the set { F r }. Here we are inter¬ 
ested in the dynamics of fluctuations of F(f), defined in the 
equations Q (the index q has been momentarily dropped 
to simplify general formulas), due to thermal hopping be¬ 
tween different states of the system in equilibrium. This dy¬ 
namics is contained in the autocorrelation function of F 
defined by the equation (|4ji. In a stationary random pro¬ 
cess (F(r)F*(f + t)) = (F(t — f)F*(r)) which implies 
<!>(—£) = <!>*(£) therefore it will be convenient to restrict the 
analysis to the case t > 0 and write the equation ([5]) as 


The expression for l/T) involves the Fourier transform of 
$(£) (Eqfl2]l. This Fourier transform can be obtained by de- 
velopping G(f) in the basis of the eigenvectors of A. Let 
{M m } be such a complete set of eigenvectors and Mi m a 
matrix having these eigenvectors as columns, with the cor¬ 
responding set of eigenvalues {7™,}. The inverse eigenval¬ 
ues r m = 1 /7m are the characteristic correlation times of the 
system 31 . 

Then the equation ( |T3} can be written 

Gij(t) = (16) 


J q (uj) = / ($,(£) + c.c.)e lLJt dt (12) 

J 0 

Let us denote as P(*, t\i', t') the conditional probability that 
the system will be found in the state i at time t if we know 
that it was in a state i! at a previous time t!. For a stationary 
system this probability only depends on the time difference 
t — t', we therefore define a function (J,,/ (t) as the conditional 
probability that the system is in state i at time t if it was in i' 
at t = 0: 


The matrix G (t) = { Gikif )} is the discrete-space equivalent 
of the Green function for the diffusion equation. It describes 
both (a) the relaxation and (b) the thermal fluctuations at equi¬ 
librium. We will use (a) to find the relation between the ma¬ 
trices G(t) and A, and then (b) to calculate the correlation 
functions. 

From the definition of G(f) the relaxation of the system 
from an initial state described by populations n,(0) at time 
t = 0 towards equilibrium is 


riiit) =y^Gij(t)nj(0) 

j 


therefore, comparing this with Eq.(10i one immediately sees 
that 


G(£) = e -At (13) 

The correlation function defined by Eq.(|4| describes the 
thermal fluctuations at equilibrium. For stationary ergodic 
systems the time averaging can there be replaced by ensemble 
averaging at some fixed time e.g. r = 0. In other words we 
can write: 

(14) 

where P(i, 0; j, t) stands for the probability of finding the sys¬ 
tem in the state i at time t = 0 and in the state j at some later 
time t. Using P(A; B) = P(B)P{A\B) this probability is 
rii{t = 0 )Gij(t). Here n^t = 0) is simply the equilibrium 
population rf. Therefore the correlation function can be ex¬ 
pressed in function of G(f) as 

m = J2 FiF J n °i G v^ (15 > 

hj 


and thus the correlation function ( fl5| ) becomes 

m = FtFj n i e~ t l Tm (17) 


Calculating the Fourier transform (12 1 with (17) we obtain 
finally 

•M") = X 

i,j,m m 

J 2 (co) = ]T n°M im M-^.(J7 (2) Ff+ c.c.) Tm 


1 + 4 w 2 r^ 
(18) 


The nuclear relaxation rate 1/Ti is given by the Eqs.([6]) 
with ( p~8] > (where for practical calculations the zero eigenvalue 
i.e. T m = 00 should be excluded from the sum). Note that the 
contributions of different states i are weighted by their thermal 
occupation numbers nf. This means that high energy states 
which are not often visited do not contribute significantly. 

From Eq.(p~8]> it may seem that the number of distinct con¬ 
tributions that could a priori be distinguished in 1/Ti vs 
temperature is related to the the number of different corre¬ 
lation times T m . However, since this expression involves 
sums the largest contribution will come from 

the eigenvectors M m having the same symmetry as the sets 
{Fj q> }. The functions F ^ in are related to spherical har¬ 
monics of order 2 therefore we may expect that high-order 
diffusion modes do not contribute much to 1/Ti- For exam¬ 
ple, for the model of Eq. GD it is found that all but one eigen¬ 
vector are orthogonal to { F -^' 1 } so that the spectral densities 
© are determined by a single correlation time. This is some¬ 
what similar to the situation found for the case of an isotropic 
rotation (diffusion on a sphere ) 24 : in the expansion of the 
Green function in spherical harmonics only the terms of order 
2 contribute to 1/Ti which is then characterized by a single 
correlation time. This discussion explains thus the relative 
robustness of the standard Kubo-Tomita analysis even when 
applied to systems having several degrees of freedom. Exam¬ 
ples where such simplification does not occur will be given in 
the next section. 

The equations ([6]) and([T8| have been derived for a single 
pair of spins. For abundant spins such as protons in organic 
materials where many spins can be coupled, the approxima¬ 
tion considering only pairwise interactions is often justified by 
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the fact that 1/Ti decreases as the 6 th power of the distance 
between spins so that it is dominated by the spin pairs that are 
the closest. On the other hand, long range dipolar coupling 
is often enough to ensure a common spin temperature and a 
single-exponential relaxation. This is for example the case of 
several compounds incorporating bicyclo[2.2.2]octane rotors 
where the closest proton pairs that dominate the relaxation are 
those of the CH 2 groups forming the rotor blades 7 — 17 23 . In 
such conditions we can define an average relaxation rate i.e. 
the Eq.(|6|» summed over proton pairs in a molecule. Likewise, 
summing the products F^F^* in Eq.( 18 1 over the closest 
proton pairs of a rotor yields the contribution of a single rotor 
to the average relaxation rate. 


III. TOY MODELS 

In this section we present two relatively simple situations 
where at least two distinct correlation times are contributing 
to l/Xi. In both cases, although the variation of 1/Ti with 
temperature could also be modeled by a sum of two KT fits 
with arbitrary parameters, our approach goes beyond the KT 
analysis because it tells exactly how the model parameters are 
related to the microscopic parameters such as the heights of 
the barriers and the positions of the wells 32 . 


A. Model 1: single rotor with two barriers 



<f> 



We consider a single rotor with N wells of the same energy 
separated by alternating small and big barriers such as shown 
in Fig J2] (top) for N = 6. For numerical calculation we have 
taken a sinusodal potential with arbitrarily chosen values of 
1500 K and 3000 K for the barriers (since we are interested 
in the temperature variation it is natural to express all ener¬ 
gies in Kelvin, 1 kcal/mol~ 503 K), 10 13 s _1 for the attempt 
frequency and 50 MHz for the NMR frequency. 

Fig(2] (bottom) shows the simulations for three values of 
N = 4, 6 , 8. The two peaks, around 140 K and 280 K, cor¬ 
respond to the activation of motion through the small and big 
barrier respectively. Since the low temperature process is a 
restricted rotation the corresponding relative amplitude de¬ 
creases with increasing N (e.g. forN=8the~ 45 degree jumps 
are much less efficient in nuclear relaxation than the 90 degree 
jumps). The amplitude ratio is in general quite sensitive to the 
position of the potential minima which depend on the degree 
of anharmonic terms in the potential. 

B. Model 2: Two rotor system and correlated motion 

This model is inspired by recent works 17 23 where the rotors 
are bicyclo[2.2.2]octane (BCO) functional units embedded in 
crystalline lattice 7 10 . The nuclear relaxation generated by the 
rotating BCO rotors comes from the modulation of the dipolar 
coupling between the proton pairs of the CH 2 groups forming 
the rotor blades. The interesting situation arises when pairs of 
such rotors are close enough such that the interaction between 


FIG. 2. top : Potential for model 1 with N = 6. bottom : Temperature 
variation of 1/Ti for N = 4,6, 8. 


them adds a significant contribution to the potential barriers 
favoring a correlated motion. A system with pairs of crys¬ 
tallographic ally equivalent, interacting BCO rotors has been 
synthetized by Lemouchi et al. 17 . Proton relaxation rate study 
has shown the existence of two different activated processes of 
very different energies and it was conjectured that these two 
processes could be attributed to disrotatory and conrotatory 
modes of rotor pair rotation 17 23 . 

In the Appendix [A] we address the question whether the ex¬ 
istence of a correlated mode of motion can lead to appearance 
of a distinguished activated process seen in NMR relaxation. 
Intuitively this could be expected if we consider the relaxation 
by inter-rotor dipolar coupling, however it is less obvious if 
we only consider proton pairs sitting on the same rotor. As 
far as the considered compound is concerned the intramolecu¬ 
lar H-H distances are always the shortest 1 ' 7 therefore we only 
consider the intra-rotor spin pairJ^U. 

For a pair of interacting rotors we consider the set of 
the three equilibrium configurations as depicted in Fig(3] 
(top), with majority-majority (a), majority-minority ( b ), and 
minoroty-minority (Fig. c) occupations of the rotor positions 
in a pair determined in the room-temperature crystal structure 
as discussed by Lemouchi et al} 1 . The resulting conforma¬ 
tional map is shown in Fig |3] (bottom) where the arrows cor¬ 
respond to the transitions considered in our model. The corre- 
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FIG. 3. top : Possible configurations (a, b, c) of a pair of interacting 
BCO rotors, after Lemouchi et at.® bottom : The conformational 
map of the two-rotor system is defined by two angles (0 1 , (j> 2 ). The 
map is periodic by 120 degree due to the threefold symmetry of the 
rotors. In the discrete model used here we consider the equilibrium 
states a, b and c and transitions between them: full arrows represent 
the single rotor jumps and dashed arrows the correlated jumps. 



Temperature (K) 

FIG. 4. Symbols: 1/Tj data at 50MHz (crosses) and 210MHz (tri¬ 
angles) from Lemouchi et al. 17 . Solid lines: simulations using the 
model of Fig[3](see text). 


value of A„/,, that is of the process of breaking the major¬ 
ity configurations a. This assignment seems to be supported 
by Carr-Parinello simulations® showing a relatively long res¬ 
idence time in the a configurations compared to that of other 
states. As for the low temperature process characterized by 
Abb, it requires quite a low value of the energy of b configura¬ 
tion (Eb — E a ) to account for the observed 1/T1 data. Such 
a low value is however not consistent with the estimation of 
the occupation rate of b sites inferred from X-ray studies 17 . 
At any rate this simple model, where only rotor-rotor interac¬ 
tions play a role and the crystalline potential is neglected, is 
probably oversimplified so that the simulations shown in Fig|4] 
should be considered as toy model simulations rather than a fit 
to the data. 


lated jumps are those between neighboring b states where both 
rotors exhibit a gear-like rotation by 60 degrees. Here we are 
studying the simplest model where the potential barriers are 
only due to rotor-rotor interaction, this model is parametrized 
by three energy barriers A a b, Abb, A bc (the map would have 
lower symmetry in the general case where the crystalline en¬ 
vironment of each rotor is taken into account). 

In Fig|4]we show the results of simulations. The solid line 
shows the best fit with the experimental data of Lemouchi 
et al. 17 obtained with Eb — E a = 140K (0.28kcal/mol), 
E c — E a = 1180K (2.35kcal/mol), A ab = 3170K (-6.3 
kcal/mol), A bb = 1000K(-1 kcal/mol), A bc = 1500K(-3 
kcal/mol) and the attempt frequencies of 1.9 10 12 s -1 and 
8.2 10 11 s _1 for the uncorrelated and correlated jumps respec¬ 
tively. We could not find any significantly different set of 
parameters reproducing this data. Within this set the high 
temperature activated process would be related to the high 


IV. CONCLUSIONS 

In this work we have proposed a simple and general ap¬ 
proach to calculate correlation functions involved in nuclear 
spin-lattice relaxation for a system undergoing thermal hop¬ 
ping between a number of potential wells. We argue that 
there should be many cases where this approach can yield a 
much more complete information about the dynamics of mo¬ 
tion compared to that obtained from a simple Kubo-Tomita 
fit. Examples of this approach in the field of molecular rotors 
are studied where the analysis of the temperature dependence 
of the relaxation rate may provide detailed information about 
the form of the conformational potential. In particular, the 
method is capable of modelling the effects of gear-like motion 
in coupled rotor systems. The use of this approach might be 
generalized if numerical tools are developed similar to pack¬ 
ages widely used to study dynamics by simulation of NMR 
spectra. 
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Appendix A: Relaxation by intramolecular spin pairs and 
correlated rotors 


The Tj due to the intramolecular spin-spin interaction ob¬ 
viously probes only the single rotor correlation function it is 
therefore interesting to ask whether it can reveal the existence 
of a correlated mode of rotation. It will be convenient to de¬ 
fine a composite index (*, it) where i and it number the posi¬ 
tion the left and right rotor respectively. Then, considering the 
relaxation of a proton pair on the left rotor, the involved ge¬ 
ometrical factors are independent of it: Fj eit = which 

implies that the correlation function will involve the effective 
G matrix for the left rotor, defined as a contraction of G over 
the indices of the right rotor: 


/-deft 

U ij 


i',j' 


Likewise we can write the master equation for the effective 
population of the left rotor: 

dnf ft d ^ ^ 

Then, if we consider a model where the hopping probability 
of the left rotor does not depend on the current position of the 
right rotor, that is if we only consider hopping between crys¬ 
tallographic ally equivalent states (e.g., in our case 120 degree 
jumps) then AmOCw') does not depend on j'. In such 
case, whether correlated jumps exist or not, the effective mas¬ 
ter equation describes a single rotor problem with an effective 
hopping rate: 


dn x f 

dt 



7 , 


nf\t) 


This means that, even if a correlated mode of rotation ex¬ 
ists, it does not necessarily lead to distinct activated process 
that could be seen in the temperature dependence of 1/Ti. In 
our case, such correlated mode of rotation between the neig- 
boring a states would just renormalize the single-rotor cor¬ 
relation time without introducing any new correlation times. 
Therefore, if the model is to account for a distinct and observ¬ 
able process related to the correlated motion we must take into 
account all non-equivalent configurations so that the hopping 


probability of a rotor depends on the current position of the 
neighbouring rotor, as in the set of states shown in Fig [3] 
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